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SYMBOLS 


English  Nomenclature 

a,  b 
c,  c„ 

^VD 

D„D 

E,,  Ey,  E,,  E 

E'=Ey/E, 

Gy,,  G^,  G 
h 

K„,K,,Ky,K^ 

L=T-U 

M 

nip  =  Pph 

9  9 

n;,  n; 

a, 

T 

t 

U 

C/cc 

u,  V 
w 

jc,  y,  z 

Zp(x,y) 


plate  streamwise  and  spanwise  lengths  (or  length  and  width) 

speed  of  sound  in  disturbed  and  undistinbed  flow 

generalized  coordinates  in  x,  y,  and  z  directions 

viscous  damping  coefficient 

orthotropic  and  isotropic  plate  bending  stiffness 

orthotropic  plate  shear  stiffness 

orthotropic  and  isotropic  Young’s  moduli 

ratio  of  in-plane  orthotropic  moduli 

orthotropic  and  isotropic  shear  moduli 

plate  thickness 

curvature  in  n,  x,  y,  and  xy  directions 

Lagrangian 

Mach  number 

plate  mass  per  unit  area 

slope  in  n,  x,  and  y  directions 

applied  in-plane  loads  in  x  and  y  directions 

generalized  nonconservative  forces 

kinetic  energy 

time 

elastic  energy 
freestream  velocity 

in-plane  displacements  in  x  and  y  directions 
transverse  displacement 
spacial  coordinates 
initial  plate  shape  function 


Greek  Nomenclature 

a,,  Uy,  a.,  a  orthotropic  and  isotropic  coefficients  of  thermal  expansion 

P  =  (M^  -  1 )  ^  compressibility  correction  factor 

a,  ,  Py  in-plane  modal  functions  in  x  and  y  directions 

Jxy,  Jyz,  Ixz  shcaT  stTaios 

Af,  EPs  static  temperature  and  pressure  differentials 


IV 


SUMMARY 


The  equations  of  motion  of  a  constant  thickness,  rectangular  panel  subject  to  external 
static  loads  and  unsteady  supersonic  aerodynamic  forces  is  derived.  Both  isotropic  and 
orthotropic  panels  are  examined.  The  static  loads  come  from  in-plane  forces,  a  temperature 
differential,  and  a  pressure  differential.  The  equations  allow  for  slight  initial  curvature  due  to 
manufacturing  imperfections.  The  equations  are  derived  using  a  nonlinear  strain-displacement 
relationship  which  permits  both  linear  and  nonlinear  panel  response  studies.  First  ord,  r  piston 
theory  predicts  the  unsteady  aerodynamic  forces  due  to  the  supersonic  flow  over  the  upper 
surface. 


SECTION  I 


Introduction 


The  equations  of  motion  for  a  panel  are  required  to  predict  its  response  to  various 
external  loading  conditions.  The  type  of  analysis  to  be  performed  is  dependent  upon  the 
applied  loading  conditions.  Some  of  the  possible  loads  include  in-  and  out-of-plane  force 
distributions  and  temperature  and  pressure  differentials.  The  equations  developed  here  will 
allow  for  any  combination  of  these  loads. 

The  equations  will  be  used  in  a  panel  flutter  analysis  in  which  the  panel  may  be  subject 
to  all  these  loads  simultaneously.  The  in-plane  loads  are  assumed  constant  along  a  side.  The 
temperature  of  both  the  panel  and  its  supports/boundaries  are  assumed  different  and  time 
independent.  An  applied  static  pressure  differential  exists  between  the  upper  and  lower 
surfaces  of  the  panel.  And  finally,  a  fluid  flows  over  the  upper  surface  of  the  panel  at 
supersonic  speeds  (3  <  M  <  5)  which  induces  unsteady  aerodynamic  forces  on  the  panel.  The 
aerodynamic  forces  are  superimposed  on  the  static  pressure  differential.  Structural  damping 
is  included  to  model  the  energy  absorption  characteristics  of  the  panel. 

The  equations  of  motion  are  derived  for  both  isotropic  and  orthotropic  panels.  Portions 
of  the  isotropic  equations  have  been  derived  in  other  documented  efforts.  The  isotropic 
equations  were  derived  in  Reference  12  without  the  applied  static  loads,  initial  panel  shape,  and 
structural  damping.  That  report  stemmed  from  research  performed  by  the  second  author  and 
published  in  Reference  4  which  included  the  initial  panel  shtqje  and  thermal  loads.  This  and 
future  efforts  will  extend  and  validate  their  research. 

In  Section  II,  the  orthotropic  constitutive  relations  are  presented  and  the  strain- 
displacement  relations  are  derived.  These  relations  are  used  in  the  development  of  the 
nonlinear  equations  of  "»otion  for  panels  with  .flight  initial  curvature.  The  nonlinearity  is 
introduced  in  the  relationships  between  strain  and  plate  deformation  which  are  summarized  at 
the  end  of  Section  II. 

Section  III  develops  the  equations  of  motion  for  an  orthotropic  panel  subject  to  some 
unspecified  nonconservative  forces.  The  strain-displacement  and  constitutive  relations  from 
the  previous  section  are  used  to  define  the  total  strain  and  bending  energies  in  the  panel. 
Rayleigh-Ritz  approximate  modes  are  substituted  into  the  energy  expressions  which  are  then 
nondimensionalized.  The  nondimensional  energy  expressions  are  substituted  into  a 
nondimensional  form  of  Lagrange’s  equation  which  yields  an  equation  of  motion  in  each 
coordinate  direction. 

The  isotropic  equations  of  motion  are  quickly  summarized  in  Section  IV  following  the 
orthotropic  development  in  Section  III. 

A  brief  discussion  of  some  possible  transverse  acting  nonconservative  forces  is  provided 
in  Section  V.  Here  first  order  piston  theory  aerodynamics  are  combined  with  a  static  pressure 
differential  and  a  structural  damping  model  to  form  the  nonconservative  forces  acting  on  the 
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panel.  Other  aerodynamic  theories  and  structural  damping  models  could  be  substituted  in  their 
place.  Transverse  forces,  other  than  those  specified,  also  could  be  incorporated  here. 

Section  VI  presents  a  derivation  of  piston  theory  starting  from  the  conservation  of 
momentum  and  mass  equations.  The  assumptions  made  in  deriving  the  theory  and  its 
limitations  are  stated.  First,  second,  and  third  order  piston  theory  expressions  are  derived  by 
expanding  the  large  amplitude  simple  wave  pressure  relation.  A  fourth  piston  theory 
expression  is  stated  that  expands  the  range  for  which  piston  theory  can  be  applied. 
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SECTION  II 


Stress/Strain/Displacement  Relations 

The  stress-strain  and  strain-displacement  relations  are  required  to  derive  the  equations  of 
motion  of  a  vibrating  plate.  The  four  basic  assumptions  in  this  section  are:  (1)  the  material  is 
ortliotropic  with  its  principle  axes  aligned  with  the  coordinated  system,  (2)  the  plate 
aeformations  do  not  exceed  the  proportional  limits,  (3)  sections  that  are  plane  and 
perpendicular  to  the  median  plane  prior  to  loading  remain  so  af^er  loading,  and  (4)  the  material 
properties  are  independent  of  temperature. 

The  constitutive  relationship  for  an  orthotropic  material  in  a  three  dimensional  stress  state 
[Ref  6]  is 


A  thermal  strain  develops  within  the  material  when  a  uniform  temperature  differential 
exists  between  one  point  of  the  material  and  another  (or  between  the  panel  and  its  boundaries). 
The  thermal  strain  is  proportional  to  the  coefficient  of  thermal  expansion  for  the  material  and 
the  temperature  differential, 

=  aj^T  (2) 

where  j  =  x,  y,  and  z.  These  thermal  strains  are  added  to  the  right-hand  side  of  Equation  (1). 
Note,  there  are  no  thermal  shear  strains. 

The  combination  of  Equations  ( 1 )  and  (2)  is  used  to  form  the  equations  of  motion  of  a 
heated  structure  when  its  physical  dimensions  are  of  the  same  order.  This  analysis,  however, 
examines  thin  plates,  i.e.,  plates  whose  thickness  is  much  smaller  than  its  length  and  width 
{h  <  a  and  h  <  b).  Therefore,  the  plate  can  be  assumed  to  be  in  a  state  of  plane  stress  where 
cr.,  x^,  and  are  neglected.  This  assumption  is  valid  if  the  upper  and  lower  surfaces  are 


3 


assumed  free  of  stress  [Ref  2 ]  anc  ihe  compression  strain  perpendicular  to  the  median  plane 
is  small  compared  to  the  normal  strains  parallel  to  the  median  plane  [Ref  11].  The  resulting 
two  dimensional  stress-strain  relations,  written  as  functions  of  the  strains,  are 


o_  = 


(l-vi,  £,/£,) 
1 


£e+v  E  e  -  (E  a  +  V  E  a  )^T 

XX  xy  y  y  '•  x  x  xy*^y  y-*" 


£€+v  £e-(£a+v  E  a  )diT 

y  y  zy  y  X  ^  y  y  ^y  y 


T  =  G  y 

xy  xy  •  xy 


(3.1) 

(3.2) 

(3.3) 


The  stresses  are  related  to  the  displacements  of  the  plate  through  the  strain-displacement 
and  constitutive  relations.  The  derivation  of  the  strain-displacement  relations  follows  that  of 
Vol'mir  [Ref  11]  and  Timoshenko  [Ref  10].  The  normal  and  shear  strain-displacement 
relations  are  developed  for  specific  applied  conditions  and  then  all  the  relations  are  combined 
at  the  end  of  this  section.  The  displacements  in  the  x,  y,  and  z  directions  are  defined  by  the 
functions  u(x,>’, /),  v(x,  >-./),  and  w(x,y,t)  respectively.  Note  that  the  displacement  functions 
are  dependent  on  time  and  the  location  on  the  panel. 

j  Figure  1  shows  an  element  (A,  B,  C,  D)  of  the 

midplane  of  the  plate  before  and  after  deflections  are 
applied  in  the  x-y  plane.  Table  1  lists  the  coordinates  of 
points  A,  Ai,  B,  £„  A  and  £>,  where  the  subscript  1 
_  denotes  the  comers  of  the  deformed  element. 


Table  1 .  Coordinates  of  three  comers  of  undeflected  and 
deflected  in-plane  element. 


Figure  1  Midplane  element 
undeflected  and  deflected  in  the 
X  and  y  directions. 


A(0,0) 

/4](u,v) 

B(dx,0) 

BAu  +  dx  +  —dx,  \+—dx) 
dx  dx 

DiQ,dy) 

+  ^dy ,  V  +  dy  +  ^dy) 
dy  dy 

The  normal  strain  in  the  x  direction  is  found  by  examining  side  ITB.  The  infinitesimal 
length  before  the  applied  deflection  is  dx.  The  length  afterwards  is 

1/2 


A,B, 


dx^^dxW 
dx 


dv. 
—  dx 
dx 


=  dx 


1  + 


du 

dxj 


(  V'  \  ^ 

+  I  - 

dx] 


1/2 


A  binomial  series  expansion  of  the  expanded  square  root  term  yields 
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1  ^ 

1  +  2 —  + 
dx 


3u 

dx 


dy 


1/2 


,  3u  l( dy 
^  1  +  —  +  -  — 


dx  2\dx 


If  the  cubic  and  higher  order  terms  are  neglected,  then  the  relative  strain  in  the  x  direction  is 

du 


XjB, -AB 
AB 


I  (  dy 
+  — 


dx  2 


dx 


Similarly,  the  expression  for  the  relative  strain  in  the  y  direction  is 


/ 

s  = 


dv  ^  1 
^2 


dy 


du 


(4) 


If  u  and  V  are  assumed  small  compared  to  w,  then  their  derivatives  are  also  small  and  the 
squares  of  the  derivatives  are  smaller  yet.  Thus,  the  squared  terms  are  neglected  also.  The 
remaining  expressions  define  the  relative  normal  strain  in  the  x  and  y  directions 

/  du  /  dv 

e,  =  —  €„  =  — 

dx  ^  dy 

Shear  strain  is  a  measure  of  the  deformation  of  a  cubic  element  into  an  oblique  parallele¬ 
piped.  It  is  defined  as  the  angle  through  which  a  side  of  the  element  is  rotated  when  shearing 
forces  are  applied.  The  shear  in  the  midplane  element  in  Figure  1  is  the  sum  of  the  two 

angles:  0,  the  angle  between  and  the  x  axis,  and  (j),  the  angle  between  and  the  y 
axis.  The  angles  are  written  in  terms  of  the  deflections  as 


0  =  tan 


-1 


dv 

du 

dx 

(()  =  tan'* 

dy 

,  du 
dx. 

1+^ 

[  dy\ 

Multiplying  the  numerator  and  denominator  of  0  by  ( 1  -  — )  and  (j)  by  ( 1  -  — )  and 


dv  \ 


dx 


neglecting  the  squared  partial  terms  yields 

0  =  tan’^f— 

djc 


dy 


0  and  (j)  rewritten  using  trigonometric  series  expansions  are 

3  V  dj: 


0  = 

dx 


(j)  =  tan'^f  — 
‘dy 


d,  = 

dy  3  V  dy 


Only  the  first  term  from  each  expansion  is  retained  because  the  other  terms  are  much  smaller 
in  comparison.  Thus,  the  relative  shearing  strain  due  to  u  and  v  is 

du  dv 


y  xy  ^ 


dy  dx 


(5) 


Forces  applied  in  the  z  direction  cause  the  plate  to  bend,  therefore,  bending  strain 
expressions  must  be  included.  Figure  2  shows  a  schematic  of  a  midplane  element  before  and 
after  point  C  is  deflected  in  the  z  direction.  In  this  case  the  plate  does  not  move  within  the  x-y 
plane.  The  curvature  of  the  element  in  a  general  direction  must  be  defined  in  order  to  derive 
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*  z 


y 


Figure  2  Midplane  element  in  equilibrium  and  in 
bending. 


the  bending  strain  expression.  The 
slop)e  at  relative  to  A  along  the  x  axis 
is  defined  using  the  above  trigonometric 
series  expansion. 


m 


X 


dw  _  \  l  8w 
dx  3  (  8j:  j 


(a) 


the  slope  at  Dj  relative  to  A  along  the  y 
axis  is 


8w  1  /  8w'\^ 


(b) 


and  the  slope  at  C  relative  to  A  along  an  arbitrary  axis  n  is 

8jc  By 
m„  =  m, —  +  — 

"  '8/1  ^  dn 


or 


=  m^cosx]  +m  smti 


(c) 


where  r\  is  the  angle  between  the  x  and  n  axes. 


The  local  curvature  is  defined  negative  if  the  point  of  evaluation.  A,  is  deflected  convex 
downward  as  shown  in  the  above  figure.  (A  cup  in  this  oriertaticn  will  hold  water.)  The 
general  curvature,  K„ ,  of  the  middle  surface  is  defined 

K  =  --^"1  (d) 

"  8/1 


where  by  the  chain  rule 


8  8  8  . 

—  =  — cosq  +  — smt] 

8/1  8jc  8y 

The  minus  sign  was  introduced  in  Equation  (d)  because  the  second  derivative  of  w  with  respect 
to  n  is  negative  for  a  plate  with  positive  curvature.  Substituting  the  derivative  expression  and 
(c)  into  (d)  produces 


2 

- cos'^n  + 

8;c 


dniy 

dy 


sin^q  + 


+ -  smn  COST) 

dy  ) 


And  finally,  substituting  in  the  slope  expressions,  (a)  and  (b),  gives  an  expression  for  the 
curvature  along  n 


K. 


8^w 

dx^ 


8lj 


21 


COS^T)  - 


8^w 

8y^ 


4(17)^: 


sin^T) 


8^w 

8jc8y 


J. 

4 


sinlq 


The  X,  y,  and  xy  curvatures  are  extracted  from  this  equation.  The  squared  terms  within  the 
brackets  are  neglected  because  they  are  two  orders  of  magnitude  smaller  than  the  other  terms. 
The  curvatures  are  thus  defined  as 


6 


=  - 


a^w 

ay  2 


a^w 

aacay 


(e,f,g) 


The  strain  in  the  plate  due  to  bending  is  proportional  to  the  curvature  and  the  distance 
from  the  median  plane  as 

6,.  =  CzK^  (h) 

where  C  is  1.0  when  i  is  :c  or  y  and  C  is  -2.0  when  i  is  xy.  Substituting  Equations  (e),  (f),  and 
(g)  into  (h)  yields  the  kinematic  relation  between  the  out-of-plane  deformation  and  strain 


€,  =  -z- 


a^w 


// 


a^w 


-2z 


a^w 


(6) 


ax2  '  ay^  dxdy 

The  midplane  of  the  plate  stretches  when  there  are  large  out-of-plane  deflections.'  In 
this  case,  the  new  length  of  ZS  in  Figure  2  is 


AB, 


{dx)^  +  \—dx 

'  '  ajc 


1/2 


After  factoring  out  dx  and  replacing  the  square  root  with  the  first  two  terms  of  its  binomial 
series  expansion,  AB^  is 

[i .  if 


AB^  “  dx 


2\  ax 


21 


(7.1) 


The  relative  normal  strain  in  the  midplane  in  the  x  direction  due  to  large  w  is  thus 

///  _  1  /  aw)^ 

~  2  i  ax  J 

The  expression  for  the  relative  normal  strain  in  the  y  direction  due  to  large  w  is  derived  in  the 
same  manner  to  be 


///  _  1  ( aw)^ 
"  2  [  ay 


(7.2) 


The  shear  strain  due  to  large  z  deflection  is  defined  by  examining  the  length  of  each  side 
of  the  deflected  plate.  According  to  Vol’mir 


aw 


aw 


m^HdyY^  ^dy-^dx 


1/2 


dy  '  dx 

Applying  the  cosine  law  introduces  the  included  angle 

and  replacing  the  cosine  term  with  the  first  term  of  its  series  expansion  yields 


'  Deflections  are  limited  such  that  the  curvature  is  less  than  0.12  radians  per  longest  panel 
length.  Expressions  (e),  (f),  and  (g)  must  include  the  neglected  terms  if  the  curvature  is  greater 
than  this  value. 
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(B,D,f  =  {AB,f*{AD,f-2AB,AD,y,"' 

Substituting  in  the  appropriate  length  expressions  and  neglecting  the  squared  derivatives  on  the 
left-hand  side  results  in  an  expression  for  the  relative  shear  strain  in  the  midplane  due  to  large 
out-of-plane  displacement, 

///  ^  aw  ^  .7  3) 

ajc  dy 

If  the  plate  is  not  flat  before  any  load  is  applied,  then  the  initial  shape  of  the  plate  needs 
to  be  included  in  the  strain-displacement  relations.  Zp  is  a  function  of  x  and  y  that  describes 
the  initial  shape  of  the  midplane  of  the  plate,  and  w  is  the  deflection  in  the  z  direction 
measured  from  Zp.  It  is  assumed  that  the  maximum  value  of  Zp  is  less  than  or  equal  to  the  plate 
thickness.  The  total  deflection  from  the  x-y  plane  is  written 


Wt  =  Zp^  w 


The  strains  in  the  midplane  are  derived  as  before.  The  initial  length  of  side  ITE  is 

1/2 


AB  =  dx 

Using  a  binomial  series  expansion  again  yields 


and  the  length  after  the  plate  deflection  is 


^ZpV 


dx 


[i*- 

(dz.\ 

dx 

P 

f 

2 

[axj 

AB^  «  dx 


=  dx 


2\dx 


dx  dx  2 


3w 


The  relative  strain  in  the  x  direction  due  to  a  deflection  w  above  a  non-flat  plate  is 

_  _  ^  1  1 0w\^ 

iv  _  '^^3  _  dx  dx  2  \  dx  j 

1 

Multiplying  the  numerator  and  denominator  by  1  -  neglecting  the  fourth  order 

partial  derivative  terms  yields  a  new  strain  term  in  addition  to  (7.1) 

^  ,v  ^  ^  (8.1) 

'  dx  dx 

The  relative  strain  in  the  y  direction  under  these  conditions  is 

?5p  ^  (8.2) 


dy  dy 

The  shear  strain  relation  is  again  based  on  the  change  in  the  angle  between  dx  and  dy. 
In  this  case,  the  initial  and  final  angles  are 
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(8.3) 


0.  = 


It 


0,  =  -  - 
^  2 


dx  dy 

dWj.  dwj. 

dx  dy 


The  shear  strain  term  due  to  a  deflection  w  in  the  z  direction  above  a  non-flat  plate  is 


Y  •''  =  dv/ 

^  dx  dy  dy  dx 

The  strain-displacement  relations  for  a  non-flat  plate  that  is  being  stretched  in  the  x  and 
y  directions  and  bent  in  the  z  direction  are  the  sum  of  the  previous  strain  expressions 


e. 

3u 

-  z - 

1 

+  — 

'  dw'' 

dw 

(9.1) 

X 

~  lx 

dx^ 

2 

.  dxj 

dx 

dx 

dv 

d^w 
-  z - 

1 

+  — 

'  dw'l 

2^  dZp 

dw 

(9.2) 

y 

dy 

dy^ 

2 

.  dy  ] 

dy 

dy 

dv  ^ 

du  ^ 

dw  dw 

-  2z 

d^w 

+ 

1 

dw  ^  dZp  difi 

(9.3) 

Jx 

dy 

dx  dy 

dxdy  dx 

dy  dy  dx 

9 


(This  page  is  blank.) 
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SECTION  m 


Orthotropic  Equation  Development 


The  equations  of  motion  of  an  orthotropic  panel  are  derived  using  Lagrange’s  equation, 

d 


dt 


=0 


(10) 


where  are  the  generalized  coordinates.  The  Lagrangian,  L,  is  the  difference  between  the 
total  kinetic  energy,  T,  and  the  total  elastic  energy,  U.  The  total  linear  elastic  strain  energy  in 
a  body  is  defined  as 


(/=  —  //  /  06+06+0€+T  Y  +T  Y  +T  Y 

2  I  I  I  **  ^  ^  yz’yz  «*** 


dxdydz 


(11) 


and  the  total  linear  kinetic  energy  is 

T  = 


i  f  fm  f— 


dxdy 


However,  since  the  in-plane  displacements,  u  and  v,  are  assumed  small  in  comparison  to  w, 
the  kinetic  energy  will  be  written  as 


'■ilHW 


dxdy 


(12) 


The  plane-stress  assumption,  i.e.,  the  stresses  through  the  thickness  are  assumed  small 
compared  to  the  in-plane  stresses,  was  introduced  in  the  last  section.  Incorporating  this 
assumption  into  the  elastic  energy  expression  yields 


€  +  T  Y 
y  y  ^y'^y 


dxdydz 


(13) 


Applied  in-plane  loads  (no  shearing  load)  are  included  through  the  stress-sti'ain  relations, 
(3.1)  and  (3.2),  as 


- - [e  +£*v  €  -(a  +£*v  a)^T]  +  -N‘‘ 

.  2  L  *  *>  y  ^  ^  ^y  y'  \  i,  ^ 


a  = 

y 


— - fe+v  e-(a+v  a)AT 

,*  2  L  y  *y  *  '  y  xy  ^^ 


(i-£*v:,) 

The  applied  loads  are  positive  if  the  plate  is  in  tension. 


4v 


(14.1) 

(14.2) 


Equations  (3.3)  and  (14)  are  substituted  into  (13)  to  produce 


U  = 


2(l-£’ 


el*  £ •  ej *  2E •  V  6.€, ^  ( 1  -£ ■  vj,)  y", 

Ex 


-  +  E*  Ar(x,y) 
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Note  that  the  temperature  differential  with  respect  to  the  panel  is  a  function  of  x  and  y. 
The  strain-displacement  relations.  Equations  (9),  are  rewritten  as 


^  dx^ 


where  and  are  the  stretching  strain  terms.  Substituting  (16)  into  (15)  gives 


2(l-£:*vM 


eL  -  £  ’  S,  .  2£  •  .  (1  -  £  •  v^,)  Y 


■y  G  f 

•..2  \  Jty.2 


-fe  +£*v  -^€  +£*-^6  +£*v  €  la  Ar(x,y) 

I  10  ly  yol  x"* 

V  ®x  “x  ) 


d-E’v") 


(N/ £„*<£„) 


[I  Sx^  I  I  J  "  dx‘  dy^ 


+  4(l-£ 


.  2  ^xyfa^W  V] 
£,  idxdyj  . 


_  d^w  (  d^w  d^w 

-2z  €  - +  £  €  - +  £  V  e  - +e  - 

^‘’dx^  ^“dy^  n'‘’ay^  ^"ax^ 


+  2(1-£-v")^Y 

”  E,‘’‘dxdy  2(, 

®y  a^W  ®y  d^W  p, 

+  £  - +  E  - +£ 

’'a,  dx2  a,  ay2 


2£,A  I  3:t2  ’  83.2 


•v 

'"dy^J  ^ 


dx  dy  dz 


Ar(x,y) 


h  h 

The  result  after  integrating  with  respect  to  z  from  “y  to  ^  is 


((\€l*E-e%*2E- 
2(1-£'vL)  JJ  I 


V  €  €  +(1“£*V^)  — — 

^xy^xo^yo  ^  ^  ^xy/  ^  *  xy 


I  a  a  \ 

-  6  +£*v  -^€  +£*^e  +£*v  6  a  ATfx.y) 

I  xo  xy  _  xo  yo  xy  yo  x 


~xo  xy  '■xo 

“x 
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(i-£’vb 


W\{ 


»2£-v 

"  dx^  dy^ 


where 


12(l-£-vJ,) 

The  first  double  integral  is  the  total  stretching  strain  energy,  Ug,  and  the  second  is  the  total 
bending  strain  energy,  Ug.  They  are  restated  here  with  the  stretching  strain  terms  leplaced 
with  their  respective  strain-displacement  relations 

(/, - ^ + 

[-£*vL)  J  J  ^\dx)  dx 


dy  2  ^  dy  j  dy  dy 


du  I  ( dw'i^  dwlfdv  1  ( dv/Y  5w 


'’'[djc  aidjcj  Bx 


Bx  dxj  dy  2 


UyJ 


dy  dy 


£,  [dx  dy  Bx  By  Bx  By  By  Bx 
.  Udu.  Wdw'\^.  dwl  ^  , 


Aru,) 


\  du 

£,A 

[dx 

(1-E‘vl,) 

dv 

E^h 

By 

dy  2  \  dy ;  By  By 


£“  +  -L  ^ 

dx  2 1  dx  1  dx  dx  * 


By  2  V  dy  J  dy  dy  J  J 


'  2jJl{axV  (dy^J  ^’ax^dy^ 


Note  at  this  point  that  the  applied  conservative  loads  do  not  contribute  to  the  bending 
strain  energy.  However,  they  would  if  the  panel  thickness  was  asymmetric  about  the  median 
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plane. 


A  Rayleigh- Ritz  approach  is  used  to  represent  the  three  displacements,  u,  v,  and  w,  with 
the  product  of  temporal  and  spacial  (modal)  functions 

u  =  aj^(r)a,.(x)p^.(y) 

V  =  (20) 

w  =  c„,(r)(|>„(x)i|r„(y) 

where  the  modal  functions  are  dependent  upon  the  plate  boimdary  conditions  only  and  the 
indices  are  summed  over  the  number  of  modes  included  in  the  analysis.  The  temporal 
functions  are  the  generalized  coordinates. 


Equations  (10),  (12),  (18),  (19),  and  (20)  are  sufficient  to  define  the  dimensional 
equations  of  motion  of  a  plate  that  is  subject  to  the  specified  conservative  loads  and  some 
unspecified  nonconservative  load,  Q^j.  However,  the  nondimensional  equations  of  motion  are 
desired.  The  following  nondimensional  parameters  are  substituted  into  the  above  equations. 


u  =  Au 


v  =  Av 


w=  Aw 


x  =  ax 


y  =  by 


z  =  hz 


Lagrange’s  equation  becomes 


AlIL 


the  total  kinetic  energy  is  rewritten  as 


the  total  bending  strain  energy  is  now 


u.  - 

^  2  ala j  , 


mi& 


+  £*  - 


dx^  dy^ 


+  4D  - 

j'y  h 


dxdy) 


2 

where  =  —^{1-E*  v^^y).  The  total  nondimensional  stretching  strain  energy  is 

*  a  J  J  \  [  doc  2  a\dx  j  a  dx  dx\  A  j  dy  2  A  ^  dy  j  bdydy 
+2£*v  —  —  +  i.— 

A  die  2  a  (  djc  j  a  dx  dx  dy  2  A  (  dy  j  bdydy 
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dy  a 


0u  h  5w  8w  h  3w  h 


3)’  dx  dy  b  dx  dy  b  dy  dx  J 
( 3wV  dw 


n 


g  6u  ^  1 
h  dx  2 


laiF 


dx  dx 


-£* 


a. 


+  V  — 

Jty ' 


\ a  dv  ^a( 8w\^  a  ^^p  dw 


h  dy  2  b\dy )  b  dy  dy 


a^ATUy) 
a.^^T{x,y) 


.  (»-g'v> 


g  8u  1  /  aw^^ 


h  dx  2 


awV 

ajj 


al  ai 


^  ^  1  aw^^  ^  g  3w 

ay  ay 


iV  "  I  dx  dy 


E^h  b[hdy  2b[dy 

The  nondimensional  displacement  functions  are 

u  =  g)/T)a,.(x)p^.(y) 

V  =  fc«(t)at(f)p,(50 

^  ”  ^mn 

The  kinetic  energy  written  in  terms  of  the  nondimensional  displacement  functions  is 


(25) 


(26) 


J*  =  ^  ^ 

2  g 


Dxb(hV^  ^ 


(27) 


The  bending  strain  energy  written  in  terms  of  the  displacement  functions  is 


2  a[a)  yj  dx 


dx 


ay  ^  ay  ^ 


£*v  £  — —ib.iir, — ^+£*v  -  <b. — - — ^Ur. 

dx^  '  '  dp  aJ2  ay2 


ay 

+  4D 


dxdy 


b)  dx  dx  dy  dy 

And  finally,  the  stretching  strain  energy  equation  becomes 

b  CCf-  -  da^da,  h  -  -  -  5<l>t  ^<l>« 
U=6D-U  a,., a.. — - — ip.p.  +  -c.,c  a  — iur.ilr  p. 

"  ^aJJl  *'  dx  dx^^^^  a  dx  a3F 


(28) 


4lgj 


I  — I  c.,c  c  c.. — - — - — - — -ilr.ilr  ilr  ilr 

kl  mn  op  y  g-  g-  g-  g-'i^l'¥n^p'¥j 


h  -  -  dz„  a<|)j.  aa,  (hV  -  -  ^4>, 

+  2-c..,g.. — ^  —  — iur,p.+  -  c^,c..\ — ^  —  —  ilf,i|r. 

-  *'  aJc  aJc  ax^'^^  UJ  dx  ax^'^^ 
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hV-  -  -  dz.  a<i>t  a<i>„  a<t>. 


^  kt  mn  IJ  '•’/'♦’n't'y 


dx  dx  dx  dx 


+  £*  - 


ZZ  ^P,  5P>  h-  -  Z  A.  JL 

bub,, a. a, — =• — :±  +  —  o,,  <!>**«; — zr — =- — = 

^  ay  ay  b 


1  (hY-  -  -  -  ^  ^  ^  ^ 
*7  7 


ay  ay  dy 
dy  dy  dy  dy 


h-  r  .  ai-.aP; 


b  *'  dy  ‘  '  ay  ay 


( j  \  ay  j  '  ay  ay 


'A\--  r  r  ^ 


+  1—  I  C.;C  C-.  - ^  <J>|.4)  <j)* - =■  - -  - ^ 

I  I,  I  kl  mn  tj 


ay  ay  ay 


+  2£*v  - 


aa.  ap, 


-  _  1  A-  -  -  .  .  ^'ri  ^'Tn  n 

b.,a-.a. — p+ c.,c  a..d>.d> - =- — =-p. 

*'  av  ih  *'  ""  ar  dv  dv  ^ 


aa.  ai|r,  aijf, 


aj  ay  2  57  5y  5y 


1  /, - r  54>„  ,  ,  5p. 

- Ct.c  b,. - — a.ilr.Uf  — - 

2  a*'""  aj  a7  ^7 


1  /,2 «  _  —  a<j>4 


dy 

...  ^’I'p 


>  1 .2_  7  7  7  „  7,  A  A  4r,  i:;?  ,|r, 

4  a*  *'  ""  dx  dv  dv  ■' 


dx 


dx  '  dy  dy 


h  -  ^  dz  ^a^d^.  h  -  z 
_  c.,a..  — ^A. — - — -  B.  +  —  Ct.a..  — ^  — =-  a,i|fi — 
b  *'  ay  ^*  87  ay  a  *'  ax  ax  '  '  ay 


I  h^-  -  -  dz  aA.  aAi ,  aA, 


in-—  —  '^'■p  '^^*1  ,  tb 

+ - c.,c  c.. — - - A  — Ai — ^A- 


2  afc  *' '""  ay  87  "  87  '  ay 

1  /,2_  _  _  aF  aA„^  aA, ,  aA, 

- c.,c  c.. — ^At - A- - A  — ^ 

2  ail  *'  ""  a7  *  a7  '  ay  "  ay 


/,2  _  _  azp  az  aAt^  ,  aA^] 

—  C^,  c,,.  -r^  ^  -r^  A,  A,  ' 


ab 


^ki^ij 


dx  dy  dx 


dy  \ 


+  D_ 


xy 


-  -  aa.  aa,. 

L  L  K  t 


bki  bii  — ^  — =  P/  P ,  ^  I  “  I  ^ij  “t  “i 


8p,ap. 


ax  a7 


ay  ay 


a  -  -  ^Pf  h-  -  -  ^4>*  ^A« 

+  2 - ii. , a,. a, p, -^  +  2 - c. , c  „ i>. .  A„ ^ A, ^  P, 
b  *'  ax  '  '  ay  b  *'  ""  ax  "  ax  '  ay  ^ 


h-  -  -  ^  d^m 

+  2 - c.,c  a.. At - a. - A  — - 

kl  mn  ij^k  .  Q-  H'n 


+  -  C.,c__c 


Lr  V  V 

kl  mn  op 


c..<b|, — — — ^<|). - ijr  \|r  — ~ 

'^  *  a7  a7  '  ay  "  ay 
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+  2-c.,6 .  — :§<!>* — =■ — =- P,  +  2 - c^.a-  — :r4)*a, — =- — ^ 

b  ^  dx  ^  dx  dy  ^  b  b  ’  dx  dy  By 

-  A-  r  ^Z;,  ^<l>*  5  a.  .  .ah-  -  5(j)*  ap^ 

+  2  —  c.,b.  — - - tlf.P  +  2 - c.,a  — - - a  U;, — - 

b  By  Bx  Bx  •'  b  b  '■'  3y  ‘  ^  By 

T  ^  ^l'-...  ^1'. 


^{hy-  -  -  5Zp^  d<l>„  d^lf, 

+  2  —  c.,c  c  — i-<b. - u) - Uf 

*'  Bx  *  Bx  '  By 

hy-  -  -  „  ^'1'. 


By 


>2  -  c„c  c.  — £^4)  -^4,, ^4,. 

I  kl  mn  ,j  Q-  Vj 


ay 


>‘Y-  - 


>1-1  C^,C.J 


Bx 


By  By 


(hy-  -  /5zp\2a<i)t  a<i>. 

^(fc)  ay)  Bx  aj’^'’^^ 


a  -  aa,  1  _  _  a<|).  a4)i 

la., — ^P  >-c.,c..  — —  Uf.Ur. 

h  Bx  ^  2  *'  a;  Bx 


-  az  a<i>,. 

+  c..  — - - ilr. 

Bx  Bx 


Ar(J,y) 


a. 


az  1  fl-  -  X 

—  h,,  a. — i-  + - c.,c.. d).Q>. - - 

h  '  By  2  b  *'  "  ay  ay 


a_  az  anf 
b  ay  ay 


a,Ar(x,y) 


Eh 


a  _  aa  1  _  _  a<|)t  a4>,.  _  az  a<|),. 

—  a.. - B.  +  —  c.,c.. - Ur.ijr.  +  c..  — - - Uf . 

h  Bx  2  *'  aJ  aJ  Bx  Bx  ^ 


Nl 


(l-£*v,,)a 


E^h 


az  1  a  -  -  X  X 

—  b  .a.  — -  + - c.,c  ■  (b. d>. - - 

h  ‘J  •  3-  3- 


By  2  b 


a  -  az  B^ 
+  —  c..  — -w. — - 
b  By  ‘  By 


By  By 
Ny^dxdy 


(29) 


One  other  assumption  must  be  made  before  the  three  nondimensional  equations  of  motion 
can  be  derived  from  the  above  equations.  The  only  forces  acting  on  the  plate  that  have  not 
been  explicitly  defined  to  this  point  are  the  nonconservative  forces.  It  is  assumed  that  all  the 
nonconservative  forces  act  in  the  z  direction  only.  Thus,  the  in-plane  nonconservative  forces 
are  neglected. 


The  equation  of  motion  in  each  direction  is  determined  by  substituting  each  generalized 


17 


coordinate  into  the  nondimensional  form  of  Lagrange’s  equation,  (22), 


dU, 

0 

da,j 

dU^ 

=  0 


(30.1) 


(30.2) 


dl  dT] 

dt. 


+  <?  =0 
^zmn 


(30.3) 


X  Equation 

The  X  equation  is  formed  by  substituting  the  stretching  strain  energy  expression  into 
Equation  (30.1).  Note  the  kinetic  and  bending  strain  energy  expressions  do  not  have  any  terms 
containing  a.j .  The  derivative  of  Ug  with  respect  to  a.^  is 


ffh‘ 


_  da.  da.  h  -  -  ^4>*  da.  u  -  dz„  d^.  da, 

<■„  ^ ^  P,Py  *  -  ^ ^ Pj  •  2^  J.,  J ^  ^ i|r,P, 

dx  dx  ^  a  dx  dx  dx  ^  a  dx  dx  dx  ■' 


da,  d^i  i  A_  _  5a,  dilr, 

+  2£*v  -  Z»^,a. — i  — 3.  +  --c,,c  (b.d)  — 

*'  *  dx  ay  “  2  b  dx  dy  dy 


‘  ai  ay  2  b  " 
.  h- 

~Cn — — =■ ='P{ 

b  *'  ay  *  ax  dy 


a\a-  ^P/^Py  -  da^  ap.  h-  -  .  ^4>„  d}\f,  ap 

[a  ‘  ay  ay  *'  ax  ■  '  ay  b  *'  *  ax  '  ay  "  ay 

h~  dz  5i|f,  ap  h_  dz  d^^  ap 

b  ax  *  '  dy  dy  b  *'  ay  ax  'ay 


=  .  (1-eXJ 


— ^^Adxdy  =  0 
Eh  h  '  ax  ^1 


Simplifying  for  a  specific  ij  and  regrouping  yields 
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^kl  ^mn 


(aV^^nl  da^d}\f, 

Uj  3y[  *  aic  ay 


a<b.  ap. 

If  B.  +/)  _!l*a.,|,  _LZ 


+  a  Jl  +  £ 


(l-£-v,^,)^  /-aa, 


.... 


E,h  h, 

which  is  the  nondimensional  equation  of  motion  in  the  x  direction. 


^dxdyc^i 


(31) 


V  Equation 

The  y  equation  is  also  formed  by  substituting  only  the  stretching  strain  energy  expression 
into  Equation  (30.2).  The  result  is 


//( 


2£*  - 


2r- 


z  ih- -  ai|f, ai|r„  ap. 

*'  *  '  ay  ay  2  g^  g^ 

h~  ^  ai|r,  ap^.] 


+  —Cn 

b 


ay 


P  A.  ^  ‘  '  J 


dy  dy 


+  2£*v  - 


_  da.  ap.  1  _  5<l>i 5(|>  ap. 

dx  '  '  ay  2  a  *'  dx  dx  ‘  "  ay 

h-  dz  d^k 


a;c  dx  ‘  ‘ 


+  2D 


jcy 


ay  j 

-  da.  da,  a-  aa,.  ap, 

dx  dx  '  ^  b  *  aj  ay  ^ 


T  ^“'.1,  ^’^"n 

aic 


ax  ‘  ay 


b-  b-  5Zp^4)*0a;  - 

—  c., — ^<j). — =■ ^P.+  —  Cw — - — - — ^i|r,P. 

h  w  a  «  a;:  .Sv. 


ax  dx  dy 


dy  dx  dx 
2 


r*(^y^  )aa  ^  ^Py  (1  E  v  )aa  ^Py]  j-j-  ^ 

-£  — +v  -- a  Ar(x,y)  a.— ^  ^ - N^a.—±)dxdy  =  0 

Ux  nbh  '  ay  E^h  b  h  ^  ‘  dy  ] 
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Simplifying  for  a  specific  ij  and  regrouping  yields 


*kl 


'kt 


^kl  ^mn 


C  d4>t  da,  _C  dilr„  _] 

2!l({&\e:  £^a,|,3.z>  4,,£!L£*!p.l 
ftJJlaifi  ‘>  31  •^' Sy  3x  3y'^\ 

(aV  5 P,-  34>»  da.-  11  _ _ 


dy 


(l-£*vM 


which  is  the  nondimensional  equation  of  motion  in  the  y  direction. 


(32) 


z  Eoriation 

The  z  equation  is  defined  by  substituting  the  kinetic,  bending,  and  stretching  energy 
expressions  into  (30.3).  It  is  written  here  in  terms  of  a  specific  i/. 


jl^j'4>k4>idxj‘ilr,^.dyjFki 


dy^  dy‘ 


d^4>t  -f  d^^i  -  C  -Cd^'^i 

J  dx  dx  J  dy  dy 
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m 
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a<j>4  a<|)„  a<j)^ 


faV  5z„  ^z, 


b  I  dx 


X  dy  ( 


dx  dx 

a4>^  d}\fi  ^<l>t .  ,  5i|r 
ax  ay  '  dx  dy 


dxdy 


y\^ki 


( n\^  d^-  a\|f,  auf 


*2(£-v.,t2D.,) 


dx  dy  dy 


dx  ay  ^ 


ay 


aV  dZp 
b)  dy 


(aV  ^'1'-  ^4f,- 

3£*  £  4,.* 

I  L  I  a—  a~  a~ 


dy  dy  dy 

a<|)*^  a(|>..  ^  a_4^ 

ay 


^2(£’v  ^2D 


a<|).  a  A  ai|f. 


I  dxdy  c, 


t/  ^mn 


^  (f  “^S 

"Uj  J  a?  J  dy  dy'’'  ^ 
^^(bjJ  *  aJ^‘  J  ay  ay  "•" 


,C  C 
'  mn  op 


J  j  [  ax  (  ax  ax  ^  I  ^  j  *  '  5y  ay  j 


dz„(  da.  ailf.  ad).  aB,  M  _  \ 
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<>*U  UJj  J  Sy  dy 

C  ^4*,  -r^Pi  - 

/•Saj  8<|),  _ /*  0i(r, 

+  D  I — ^<b  —  dx  I  B,— ^Ur.dy 
"j  a;  r'  a?  ^ 

/*aoj  a^)-  _  /*  aiii,  _i  _ 
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auft  aiir, 

Ar(x,y)  4>*<l>,-^-^rfxrfy  c„ 
dy  dy 


r  c  dz  atb. 

ljATay)^^,^d,dy 


(l-£’v,,)/a^^ff  ra*.  «'!><. r 


E^h  i^j  Ui 


*J_£(£fo..  =  0 
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SECTION  IV 


Isotropic  Panels 


The  nondimensional  equations  of  motion  for  a  panel  made  of  an  isotropic  material  rather 
than  an  orthotropic  material  assumed  previously  are  stated  here  with  little  derivation.  The 
stress-strain  relations  for  a  thin  isotropic  panel  imder  an  applied  temperature  differential  and 
in-plane  loads  are 


Ox  = 


(l-v2) 

E 


1 


[€y+v€,-(l  +  v)aAr]  + 


(34.1) 


(34.2) 


U  = 


(l-v^)^^  ^  ^  h 

=  Gy,,  (34.3) 

Substituting  Equations  (34)  into  the  total  linear  strain  energy  expression,  (13),  yields 
E 


2(1 


b;III 


€?  +  €?  +  2ve,€y+  ^^-^Yxy-(l  +  v)aAr(x,y)(6^  +  6,) 


-X  '■> 


dxdydz 


(35) 


and  substituting  the  strain-displacement  relations,  (16),  into  (35)  gives 


eL  +  €?„  +  2  V  €„e„.  +  Yxyo  "  ( 1  +  v)  a  Ar(x,y)  (e^^  +  e^^) 


'xo  yo  ~  '  ''xo~yo 

,2 


Eh 


Uy^j  dx^  dy^  Wxdyj 

( 


+  z' 


- +  €  - 

3,2  yo  a,.2 


+  V  6. 


dx^ 


+  (1-v)y 


d^w  (1+v) 


xyo 


dxdy 


aATix,y) 


d^w  d^w 
-  +  - 


2Eh 


dx^ 


df 


II 


dx^  dy^ 
dxdydz 


h  h 

Integrating  with  respect  to  z  from  to  —  and  separating  the  result  into  the  bending 

2  2 

and  stretching  equations  using  (9)  produces 
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2v 


dx^  dy^ 


^2(l-v) 


y 

Ujcayj  J 


dxdy 


(36) 


and 


£ft  ffffau^w^) 

2 

dv  ^  1 

'^dw^  ^  ^Zp 

^  2(1-v2)  j  J  2(dxJ 

dx  dx 

dy  2 

,  dyj  dy  dy 

dx  2\dx  j  dx  dx 


dv  ^  1  I  dwj^ 


dvf 


,  (1-v) 

2 

-a(l  +  v) 
(l-v^ 

+  J: - 

Eh 


dv  du  dw  dw 
—  +  —  + - + 

dx  dy  dx  dy 

dx  dy  2  ( dx  I  2  [  dv  dx  dx 


dy  2  (  dy  j  dy  dy 

^  dw 
dx  dy  dy  dx 

dw 


3w 


dy  dy 


Ar(x,y) 


1[^  +  dZp  dw1/^a\ 

dx  2idxj  dx  dxj' 


d-v") 

Eh 


dv  ^  1  (  dw)^  ^  dw 


vdxdy 


(37) 


[  dy  2  (  dy )  dy  dy 

Equations  (12),  (36),  and  (37)  are  nondimensionalized  using  the  parameters  defined  in  (21)  but 
replacing  D,  with  D.  The  nondimensional  kinetic  energy  is 


the  bending  strain  energy  is 

flOTKS 


Us  = 


w^2 


aW  d^w)^ 


b)  Idy^j 


:n2v(£ 

b 


d^w  d^w 
dx^  dy^ 


+  2(l-v)  ^ 


2™  \21 


d^w 
dxdyj 


dxdy 


(39) 


and  the  stretching  strain  energy  is 

—  f  f  I  [—  +  —  —  ( — — 

^  i  i  1  [  2  a  (  dx  j  a  dx  dx 


Us  =6D 


dv  ^  1  ^ ( dw 


+  2v- 


b 

(1-v) 


du  ^  1  ft  ( dw 
dx  2  fl  V  dx 


^  ft  ^Zp  aw 
a  dx  dx 


dy 

dv 

[dy 


ft  dZp  aw? 
b  dy  dy 

^  1  ft  ( dwV  ^  ft  ^Zp 
2  ft  (  dy  j  ft  dy  dy 


w)^ 

2ftl57j 

-\2 


dv  ^  ^  h.§^i5L+  b  dw  ^  ft  ^Zp 

dx  ft  dy  ft  dx  dy  b  dx  dy  ft  dy  dx  J 
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dw 

Ixlx 


+ 


aw 
dy  dy 


Ar(x.y) 


+ 


(l-v^) 


a 3u ^  1  3w 
h  dx  2\ dx 


aw 

+  - £. - 


wV 

Xj 


K) 


Eh 

(1-v^)  a  a  dv  ^  la/  aw 
Eh  b[h~dy^2~b[dy)  b  dy  ayj 


ax  ax 

.T.\2  ^  dz.  dw 


(40) 


Isotropic  X  Equation 

The  nondimensional  x  equation  simplified  for  a  specific  ij  is 

•  -il/f  f 


'kl 


dz.  a<t>.  att; 

—  — — + 

ax  dx 


(1-v)  (a 


:i|  ♦..,^£^1 

h)  ‘  '  By  By 


. -  “i?i/S'"7  >.*■  ("•I 


Uj  dy 


(41) 


Isotropic  V  Equation 

The  nondimensional  y  equation  simplified  for  a  specific  ij  is 


2^ 


J  dx  •  J  dy  ^  2  J  ^dx  J  dy^^  ^ 


+  2 


2  J  dx 


Xl  da,  _/* 
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J  dy  dy  ay 

f  d(i>.  d<it  _/*  ^P,  _ 

J  dx  dx  ‘  J  dy 

f  d4>L  da.  _  f 

+  (l-v)  -^djc  ,|,,— lip.dy 

J  dx  "  dx  J  '  dy  ^ 

-2*  f  n^[v:^a  ,(,,^  .  <il!0  4,.^£*ip 
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^kl  ^mn 


+  a(l  +  v) 
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i^ljli.nx,y)a,^d;dy 


dy 


^dxdyc,. 


Isotropic  z  Equation 

The  nondimensional  z  equation  written  in  terms  of  a  specific  ij  is 

4  {/$  S  *  (f  ""'/S'' 
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(42) 


+  2 


fV*  £^1 
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-  a(l  +  v) 
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Ar(A:,y) 


1  ^4)*  54>i 


dx  dx 


dilr,  aty 


ay  dy  i 


dxdy  Cf^i 


£*i 

ay  ■  ay 


dxdy  1 


\2 


1 1  /♦**'  ^  ‘'^1  *  //5  > 


+  _L£('£)'o. 

6D  dUi  "' 


0 
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SECTION  V 


Nonconservative  Forces 


are  the  nonconservative  applied  transverse  loads.  In  this  analysis,  the  loading  is  due 
to  a  static  pressure  differential  between  the  plate’s  upper  and  lower  surfaces,  the  aerodynamic 
flow  over  the  upper  surface  of  the  plate,  and  structural  damping.  From  linear  theory,  0,,^  can 
be  written  as 

<?„/  ■  f  I {a/>,  *  (P-P.)  -  P^}^  dxdy  (44) 

where  is  the  static  pressure  differential,  (P  -P^)  is  the  supersonic  aerodynamic  pressure 
differential,  and  P^^  is  a  structural  damping  load.  The  static  pressure  differential  is  positive 
when  the  upper  surface  of  the  plate  is  in  tension. 


The  supersonic  aerodynamic  pressure  differential  is  written  as 

p.f/i 

P-P_  =  Z^LJL 


'M^-2 

1  dw] 

I  .  P-V-  3^, 

[dx 

dt\ 

1  P  dx 

(45) 


where  the  term  within  the  braces  comes  from  first-order  piston  theory  and  the  last  term  is  a 
static  pressure  due  to  the  plate’s  initial  out-of-plane  shape.  Equation  (45)  is  derived  in  Section 
VI. 


The  mathematical  model  used  to  represent  the  structural  damping  is  that  of  viscous 
damping 

Po^  =  (46) 

Substituting  the  nondimensional  parameters  into  (44),  (45),  and  (46),  and  combining  the 
results  yields 


=  hab\  I \AP,- 


9,U^  h 

dw 

4. 

a  1 

(  ) 

P.f^i  h  dZp 

P  a 

die 

uj 

dx 

P  a  dx 

+  Cyjjh 


D  \  1  dvt  —  — 


m^a 


dx  dc. 


and  substituting  in  the  nondimensional  dynamic  pressure  and  mass  ratio  expressions 


k  = 


D. 


P  = 


produces 
O  =  — * 

p 


p,g 


a[a)JJ\kD^h  ^  dx 


dx  dx 
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+  c 


VD 


f— 1 

Uj 


1/2 


8w  ]  5w 

—  dxdy 


p„t/,  at  ac 


The  final  step  is  to  replace  w  with  its  nondimensional  modal  expansion,  (26).  The 
nonconservative  forces  on  the  plate  written  in  terms  of  a  specific  ij  are 


(47) 
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SECTION  VI 


Unsteady  Aerodynamics 


Unsteady  supersonic  aerodynamic  forces  are  the  dynamic  driving  force  in  a  panel  flutter 
analysis.  Many  different  theories  could  be  used  to  describe  these  forces  but  in  this  analysis 
piston  theory  aerodynamics  will  be  used.  The  nonlinear  form  is  derived  from  the  conservation 
equations.  Several  assumptions  are  made  pertaining  to  the  flow:  (1)  the  isentropic  relations  and 
the  perfect  gas  equation  of  state  apply,  (2)  gravity  effects  are  negligible,  (3)  the  flow  is 
irrotational,  and  (4)  the  fluid  may  be  treated  as  a  continuum.  This  derivation  uses  information 
gathered  from  Shapiro  [Ref.  8  and  9]. 

The  conservation  of  momentum  equation  (including  the  irrotational  assumption)  is 

+  =  -  — VP  (48) 

dt  2  p 

Since  the  flow  is  irrotational,  its  velocity  can  be  represented  by  a  potential  function  as 

V  =  (49) 

Substituting  the  velocity  potential  into  the  conservation  of  momentum  produces 

+Iv(V$)2  = -IvP  (50) 

dt  2  p 

After  taking  the  dot  product  of  (50)  with  a  unit  vector  along  a  streamline,  integrate  from  an 
upstream  reference  point  where  the  flow  is  steady  to  some  point  over  the  panel 

•*08  OO 

which  yields 

+  fldP  =  -U!  (51) 

dt  2  J  p  2 

* 

where  is  the  uniform  and  steady  flow  velocity  at  the  reference  point. 


Since  the  flow  is  isentropic,  the  isentropic  pressure/density  relation 

p\i/r 


P  =  P 


if 

\  *  •/ 


can  be  substituted  into  the  integral  term  of  Equation  (51) 


J  p  p-  J 


p-^'ydP  = 


y-i  p. 


/  P  yr-u/Y 


Substituting  this  and  the  speed  of  sound  definition 


(52) 
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2  yP. 

=  - 


into  Equation  (51)  yields 

a* 


2 

dt  2 


Y-1 


p  \(Y-1)/Y 


2 


Also  note  that  the  isentropic  pressure/speed  of  sound  relation 


P  yY-i)/Y 

T 


(53) 


(54) 


could  be  introduced  at  this  time 


- +  _( V®)‘‘  + - 

dt  2  Y-1 


(j)’ 


1 


(55) 


The  flow  ahead  of  a  piston  in  a  tube  is  assumed  to  be  one-dimensional.  Piston  theory 
aerodynamics  assumes  that  a  vibrating  structure  in  a  supersonic  flow  is  analogous  to  an  infinite 
number  of  oscillating  infinitesimal  pistons  placed  side  by  side  forming  the  shape  of  the 
structure.  It  further  assumes  that  the  flow  changes  due  to  one  piston  do  not  effect  the  flow 
over  an  adjacent  piston.  Thus,  the  perPirbations  to  the  flow  over  a  point  on  a  panel  will  be 
defined  using  one-dimensional  theory.  Equation  (55)  written  in  one  spacial  dimension  is 

^  Y-1 


dt  2l,a2j 


i-l 


- 1 


=  if/2 
2 


(56) 


where  it  was  defined  earlier  that  z  is  the  vertical  direction.  Solving  Equation  (56)  for  the  speed 
of  sound  in  the  disturbed  flow  gives 

v2  . 

(57) 


=  c.  -(y  -1) 


-  +  I  -  I  2  12. 


dt  2(  dz 

Since  the  flow  over  the  panel  is  one-dimensional, 

Mx,y,t)  =  ^ 
dz 


(58) 


where  w  is  the  vertical  displacement  of  the  point  on  the  panel  being  evaluated. 


The  one-dimensional  conservation  of  mass  equation  written  in  terms  of  the  velocity 
potential  is 


dp  d^  dp  d^<b 
dt  dz  dz  dz^ 


(59) 


Differentials  of  the  one-dimensional  form  of  Equation  (51)  will  be  used  to  eliminate  density 
from  the  conservation  of  mass  equation.  Differentiating  Equation  (51)  with  respect  to  z 
produces 
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a®  1  dP  ^ 
dzdt  dz  dz^  P 


a^<I>  d<If  d^<t> 
-  + - 


and  solving  for  the  pressure  derivative  yields 

iL  = 

dz  dp  dz  ^[azar  az  dz^. 

The  isentropic  pressure/density  equation,  (52),  can  be  rewritten  in  differential  form  as 

dP  _  yP. 
de  pi 

which  is  then  substituted  into  (60)  to  get 


.Y-l 


i£  = 

P'  ( 

d^d>  ddf  d^df 

dz 

yPJ 

.  p  J 

dzdt  dz  dz^. 

Differentiating  the  one-dimens.'onal  form  of  Equation  (51)  with  respect  to  t  leaves 


d<b  a^O  1  dP  n 

dt^  dz  dzdt  p  dt 


and  solving  this  expression  for  the  pressure  derivative 

dP  dP  dp 

—  = - =  -p 


dt  dp  dt 

Equation  (63)  can  be  rewritten  using  (61) 

Ip  =  — bI 

dt  yP, 


a^^  ^  a®  a2$ 


.  p 


dt^  dz  dzdt 


d^df  a®  a^^ 

-  + - 


dt^  dz  dzdt] 

Substituting  (62)  and  (64)  into  the  conservation  of  mass  equation,  (59),  yields 


yP, 


.Y-l  - 


m 

dz  I 


dz‘ 


dz  dzdt  dt^ 


=  0 


Using  the  isentropic  relations,  (52)  and  (54),  Equation  (65)  is  rewritten  as 


- 


ddfV^ 
dz, 


d^  _,^ddf  d^df 


d^9 


dz" 


dz  dzdt  dt^ 


=  0 


(60) 


(61) 


(62) 


(63) 


(64) 


(65) 


(66) 


Equation  (66)  is  a  nonlinear,  homogeneous  partial  differential  equation  of  second  order  which 
describes  the  unsteady  motion  of  the  flow  over  an  oscillating  piston. 

The  theory  of  characteristics  can  be  used  to  solve  partial  differential  equations  (PDE)  that 
are  hyperbolic  and  have  the  form 

Adf^^  +  2B9^,  +  Cdf„  =  0  (67) 

where 
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dzdt 


zz 


d^<j> 

dz^ 


dt^ 


and  the  coefficients  A,  B,  and  C  are  functions  of  2,  t,  and  .  Equation  (66)  is  in  this 
form  and  is  hyperbolic  since  B^- AC  =  >  0.  The  hyperbolic  nature  of  the  PDE  dictates  that 

there  are  two  real  characteristics  that  satisfy  Equation  (66).  A  characteristic  describes  the 
motion  of  a  pressure  wave  through  a  fluid. 


The  flow  over  the  panel,  therefore,  consists  of  simple  waves  since  its  PDE  has  the  form 
shown  in  (67).  These  simple  waves  are  pressure  waves  emanating  from  the  oscillating  piston. 
Their  amplitude  may  be  small  or  large. 


The  equation  describing  the  motion  of  each  wave  must  satisfy  the  governing  PDE.  There 
are  two  families  of  waves  that  satisfy  the  PDE.  There  are  upward-  and  downward-traveling 
pressure  waves  (traveling  with  respect  to  a  fluid  particle).  The  nomenclature  assimies  the 
piston  is  oscillating  in  a  vertical,  constant  cross-sectional  area  duct  oriented  such  that  the 
upward  direction  is  positive.  Either  one  of  these  pressure  wave  families  may  exist  in  a  flow, 
but  not  both  at  one  time. 


A  flow  that  contains  simple  waves  is  defined  such  that  and  are  unique  functions 
of  each  other  [Ref.  9].  Thus,  the  derivatives  in  Equation  (67)  can  be  rewritten  as 

do,  dO, 

O  =  - 5  ^ i 

dz  do,  dz 


O 


tt 


dO, 

“dT 


do^  do, 

~df  dz 
do,  dO^  do,  dO, 

dO^  dr  dOj  dz 


and  substituting  these  into  (67)  gives 


dO, 

Solving  for  - results  in 

dO^ 


dO, 

A+2B - -^C 

d<b. 


dO.\ 


dO. 


=  0 


t) 


^  ^  -b±\Ib^-ac 

dO^  C 


and  after  substituting  in  the  values  of  the  coefficients  from  Equation  (66) 


dO, 


-(w  ±  c) 


(68) 


where  the  minus  sign  corresponds  to  upward-traveling  waves  and  the  plus  sign  to  downward¬ 
traveling  waves.  Now,  differentiating  Equation  (57)  with  respect  to  w  yields 
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(69) 


dc 

_  Y-1 

_  Y-1 

—  = 

- +  w 

- +  w 

5w 

2c 

dw 

2c 

d<»^ 

and  substituting  (68)  into  (69)  produces 

dc  Y  - 1 
dw  *  2 

Integrating  (70)  from  the  upstream  reference  point  to  some  point  over  the  panel  gives 

Y  “  1  /  •  •  V 

C  -  C.  =  ±  ■■  ^  I  w  - 


(70) 


Dividing  through  by  c«  and  noting  that  w,  =  0  because  the  flow  is  uniform  in  the  x  direction 
at  the  reference  point  yields 

A  =  (71) 

2  c. 

Recall  the  isentropic  relation  given  in  Equation  (54).  Replacing  the  speed  of  sound  ratio  with 
the  pressure  ratio  provides 

12. 

(72) 

Equation  (72)  is  written  for  upward-traveling  waves  only  since  both  upward-  and  downward- 
traveling  waves  cannot  exist  at  the  same  time.  Equation  (72)  is  the  equation  upon  which  piston 
theory  is  based  [Ref.  1  and  7].  It  relates  the  pressure  on  the  face  of  a  piston  to  its  motion 
regardless  of  the  magnitude  of  that  motion. 
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Equation  (72)  can  be  expanded  using  a  binomial  series  expansion  to  get  different  piston 
theory  orders.  First  order  piston  theory  is  generally  written  as 

=  P.c-w  (73) 


The  second  order  binomial  expansion  gives 


P-P„  = 


(74) 


And  third  order  piston  theory  is 


P-P,  = 


w  ^  y  +  1 
c.  4 


2 


+ 
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(75) 


The  expression  defining  the  velocity  of  the  piston  surface,  (58),  needs  to  be  discussed 
further.  As  stated  in  the  Symbols  section,  Zp(x,y)  defines  the  initial  shape  of  the  plate  and  w 
is  the  vertical  displacement  measured  relative  to  this  initial  sh^.  Since  the  plate  is  being 
divided  into  infinitesimal  pistons,  this  function  also  describes  the  initial  orientation  of  the 
piston  surface.  The  vertical  velocity  of  the  piston  in  a  uniform  horizontal  stream  is  then 
described  by 
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w 


(76) 


D(w*zJ  aw  dv/ 

_ £1  =  ^  fj  +  [J  _ £ 

Di  dt  “  dx  "  dx 

Substituting  the  vertical  displacement  expression  into  the  first  order  piston  theory  expression 
yields 

P-P  -  ^  aw  dw  ^  ^Zp 

A/,  if/,  dt  dx  dx 

First  order  piston  theory  is  used  in  this  analysis  since  higher  order  aerodynamics  do  not  affect 
the  panel  response  unless  the  deflections  are  on  the  order  of  the  plate  length  [Ref  5]  which  is 
outside  the  limits  of  the  structural  assumptions  imposed  here. 

Equation  (77)  can  be  employed  when  the  Mach  number  is  sufficiently  high,  M>2,  and 
the  length-to- width  ratio  is  low,  a/6  <  p  [Ref  3].  Equation  (77)  can  be  modified  to  include 
Mach  numbers  down  to  \/2. 

P-P  (U) 

P  U.  dt  dx  dx) 

This  form  was  derived  from  a  first  order  expansion  in  reduced  frequency  of  the  exact  two- 
dimensional,  unsteady,  linearized,  potential  flow  equation.  Equation  (78)  can  be  used  to  define 
the  unsteady  aerodynamic  loads  imposed  on  a  structure  due  to  a  flow  whose  Mach  number  is 
greater  than  \/2. 
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SECTION  vn 


Conclusions 


The  nondimensional  equations  of  motion  were  derived  for  a  thin  panel  using  nonlinear 
strain-displacement  and  linear  stress-strain  expressions.  The  panel  was  assumed  to  have  a 
constant  thickness  and  a  flat  or  slightly  curved  surface.  Both  orthotropic  and  isotropic 
materials  were  examined.  Solving  these  equations  would  produce  the  in-  and  out-of-plane 
motion  of  the  panel  subject  to  external  loads.  The  out-of-plane  displacements  were  limited 
such  that  the  local  ciurature  anywhere  on  the  panel  would  not  be  greater  than  0.12  radians  per 
panel  length.  Also,  the  magnitude  of  the  out-of-plane  displacements  was  assumed  to  be  much 
larger  than  the  magnitude  of  the  in-plane  displacements. 

The  boundary  conditions  on  the  panel  were  not  defined  other  than  a  fluid  flow  only  over 
the  upper  surface.  Thus,  these  equations  are  eq}plicable  to  rectangular  panels  whose  edges  are 
any  combination  of  the  following  boundary  conditions;  clamped,  simply  supported,  and/or  free. 
The  boundary  conditions  are  introduced  through  the  nondimensional  modal  functions.  Equation 
(26). 


This  derivation  assumed  that  the  panel  was  subject  to  static  and  dynamic  external  loads. 
Several  different  loading  conditions  were  examined.  The  static  loads  consisted  of  in-plane 
forces,  a  temperature  differential,  and  a  pressure  differential.  An  unsteady  supersonic  flow 
over  tlie  upper  surface  of  the  panel  was  modelled  using  first  order  piston  theory.  Structural 
damping  was  incorporated  using  a  viscous  damping  model  which  only  affected  the  out-of-plane 
motion.  The  equations  were  written  such  that  other  aerodynamic  and  damping  models  could 
be  readily  substituted  for  the  current  ones. 

The  equations  describing  the  motion  of  an  orthotropic  panel  subject  to  the  above  in-plane 
loads  and  some  arbitrary  out-of-plane  loads  were  defined  in  Equations  (31),  (32),  and  (33). 
Their  derivation  and  the  assumptions  made  during  the  derivation  were  discussed. 

Equations  (41),  (42),  and  (43)  described  the  isotropic  panel  motion  in  the  x,  y,  and  z 
directions,  respectively,  and  Equation  (47)  defined  some  possible  nonconservative  transverse 
forces  acting  on  the  panel  in  the  z  direction. 

A  large  amplitude  simple  pressure  wave  expression.  Equation  (72),  was  derived  from  the 
conservation  equatons  using  the  method  of  characteristics.  This  expression  described  the 
pressure  on  the  surface  of  a  piston  oscillating  in  a  one  dimensional  duct.  This  was  the  basis 
from  which  first  order  piston  theory  aerodynamics.  Equation  (73),  was  derived. 

This  analysis  did  not  consider  panels  with  large  initial  curvature  (shells)  or  with  an 
underlying  cavity.  Nor  did  this  derivation  discuss  the  ply  arrangement  for  a  laminated  material 
with  orthotropic  properties.  These  parameters  and  other  aerodynamic  theories  will  be  included 
in  future  research. 
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